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Abstract 

This  paper  attempts  to  relate  robust  control  and  behavioral  frameworks  by  incor¬ 
porating  structured  uncertainty  into  the  description  of  behavioral  systems.  Behavioral 
equations  are  expressed  as  linear  fractional  transformations  (LFTs)  on  an  uncertainty 
structure,  and  a  method  of  interconnection  is  outlined.  A  method  for  obtaining  input- 
output  maps  from  LFT  representations  of  behavioral  systems  is  also  presented.  This 
extension  of  the  behavioral  framework  is  compatible  with  existing  robust  control  meth¬ 
ods,  such  as  fj,  analysis,  which  can  be  used  to  provide  robustness  tests  in  behaviors.  A 
simple  example  is  presented  that  illustrates  some  of  the  issues  which  arise  in  this  extension. 


1  Introduction 

A  major  theme  in  robust  control  has  been  to  supply  the  engineer  with  a  theoretical  and 
computational  framework  that  handles  a  rich  variety  of  modeling  uncertainty  so  that  physi¬ 
cally  motivated  uncertainty  descriptions  can  be  treated  in  a  natural  manner.  In  particular, 
it  has  been  important  to  provide  computational  tools  that  analyze  systems  with  mixtures  of 
unstructured  uncertainties  and  possibly  large  numbers  of  uncertain  real  parameters.  Behav¬ 
ioral  models  are  in  turn  very  natural  when  modeling  physical  systems  from  first  principles,  or 
when  a  large  system  is  built  up  from  subsystem  models.  While  the  final  interconnected  model 
used  in  a  robust  control  design  may  have  well-defined  inputs  and  outputs,  it  is  almost  always 
the  case  that  components  are  modeled  in  terms  of  mass,  momentum,  or  energy  balances  or 
physical  laws  such  as  Newton’s  second  law,  Ohm’s  law,  and  so  on.  These  components  do  not 
have  a  pre-specified  signal  direction,  and  forcing  them  to  be  input-output  operators  assumes 
a  knowledge  of  the  ultimate  function  of  the  device.  The  reader  is  referred  to  [1]  for  a  detailed 
exposition  of  the  issues  outlined  above. 

The  behavioral  setting  provides  a  convenient  framework  for  connecting  modeling  and 
robust  control.  Uncertainty  can  be  incorporated  in  the  description  of  behavioral  systems 
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via  linear  fractional  transformations  (LFTs).  By  describing  these  systems  as  LFTs,  existing 
robust  control  methods,  such  as  //  analysis  and  synthesis,  can  be  used  to  provide  robustness 
tests  in  behaviors. 


2  Background 

2.1  Behavioral  Systems 

The  behavioral  framework  has  been  extensively  described  and  developed  by  Willems  in  [2],  of 
which  a  summary  relevant  to  this  paper  may  be  found  in  [1],  In  the  interest  of  completeness, 
a  brief  summary  follows.  A  dynamical  system  is  viewed  as  a  family  of  trajectories,  and  is 
defined  as  a  triple  S  =  (T,W,#),  with  T  a  subset  of  R  called  the  time  axis ,  W  the  signal 
space ,  and  B  a  subset  of  called  the  behavior. 

In  the  absence  of  uncertainty,  we  will  be  concerned  with  continuous-time,  linear,  time 
invariant,  finite  dimensional  dynamical  systems.  In  particular,  we  require  that  T  =  R+  and 
W  =  R?.  In  [3],  it  is  shown  that  this  is  equivalent  to  requiring  B  to  be  the  kernel  of  a 
polynomial  operator  in  the  differentiation  operator,  Jj.  We  will  also  restrict  ourselves  to  time 
trajectories  w  that  are  locally  square  integrable: 

f  ||  w  || 2dt  <  oo  (1) 

Jn 

for  all  compact  sets  C  R+,  which  will  be  denoted  w  £  Ll^ .  We  will  refer  to  the  class  of 
dynamical  systems  with  the  above  properties  as  C. 

Although  continuous  time  systems  are  considered  in  this  paper,  most  of  the  results  pre¬ 
sented  can  either  be  applied  directly  to  discrete  time  systems,  or  can  be  trivially  extended  to 
the  discrete  time  case.  The  same  can  be  said  for  the  signal  space  chosen.  We  have  chosen  to 
consider  time  trajectories  that  are  locally  square  integrable,  but  the  results  generalize  to  all 
Xp0C  spaces. 

2.2  Linear  Fractional  Transformations 

LFTs  are  a  class  of  general  linear  feedback  loops,  and  can  be  pictorially  depicted  as  follows: 


L±±_J 

X 

z 

V 

M 

W 

Figure  1:  Linear  Fractional  Transformation 


We  will  be  interested  in  cases  where  M  is  a  constant,  finite  dimensional  matrix  and  A  is, 
as  of  yet,  an  unspecified  operator.  If  we  partition  M  as 
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M  = 


(2) 


A  B 
C  D 

we  have  the  following  loop  equations: 

^  =  Ax  +  Bw  (3) 

v  =  Cx  +  Dw 
x  —  A  z 


If  the  inverse  of 
v  =  <S(A,  M)w,  where 


(I  —  A  A)  exists,  the  vectors 


and 


w 


satisfy 


S(A,M)  =  D  +  CA(I-AA)~1B  (4) 

<S(A,M)  is  called  an  LFT  between  A  and  M.  A  will  usually  belong  to  an  allowable  set  of 
operators,  A.  We  will  consider  A  structures  of  the  form 

A  =  {diag[JI,  Au]}  (5) 


where  Au  is  a  structured  set  of  bounded,  linear  operators  on  L1™,  referred  to  as  the  uncertainty 
block.  Often  we  will  need  norm  bounded  subsets  of  Au,  and  we  introduce  the  notation 

Ba,£  =  {A  G  A  :  ||  Au  ||<  e}  (6) 

Note  the  norm  bound  is  only  on  the  uncertainty  portion  of  A,  Au. 

Many  systems  with  uncertainty  can  be  represented  using  LFTs.  For  example,  uncertainty 
in  the  physical  quantities  of  a  system,  such  as  mass  or  capacitance,  can  easily  be  written  as 
LFTs  on  an  uncertainty  block.  The  same  is  true  for  arbitrary  operators  on  variables,  which 
might  arise,  for  example,  when  modeling  the  effects  of  friction  on  a  system.  The  interested 
reader  is  referred  to  [1]  and  [5]  for  examples  of  the  type  of  uncertainty  that  may  be  represented 
by  LFTs. 


3  Representation 

In  order  to  include  uncertainty  into  our  system  description,  it  is  desirable  to  first  explore  how 
systems  without  uncertainty  may  be  modeled.  Incorporating  uncertainty  into  our  system 
description  should  then  ideally  be  a  natural  extension  of  the  existing  framework. 

Consider  the  following  equations: 

x  =  Ax  +  Bw  (7) 

0  =  Cx  +  Dw 

where  A,B,C,D  are  constant,  finite  dimensional  matrices.  Then  it  can  be  shown  (see  [4]) 
that  the  behavior  of  a  dynamical  system  X  in  £  can  be  characterized  as 

B  =  {w  E  Ll2°  |  3®  €  L1™  for  which  (7)  is  satisfied}  (8) 

for  some  A,B,C ,  and  D.  Equations  (7)  are  referred  to  as  an  output  nulling  (ON)  repre¬ 
sentation  of  E,  and  is  denoted  {A,  B,C,  D}.  There  are  many  other  equivalent  state  space 
representations  of  dynamical  systems  in  £,  see  [6]  and  [7]  for  example.  As  will  be  seen, 
however,  ON  representations  can  be  naturally  extended  to  incorporate  uncertainty  into  our 
system  description. 

The  following  definitions  apply  to  the  set  {A,B,C,D}  : 
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Definition  1  ON  Representation  {A,  B ,  C ,  D }  of  E  is  said  to  be 


Dependent 

Regular 

Singular 


if,  Vs  £  C,  D  +  C(sl  —  A)-1  B  is  not  surjective 

if  D  is  surjective 

if  not  dependent  and  not  regular 


It  follows  that  regular  representations  can  never  be  dependent,  and  vice-versa.  In  the  singular 
case,  the  rank  of  D  +  C(sl  —  A)~l  B  is  normally  full  but  drops  at  s  =  oo.  It  should  be  stressed 
that  the  above  are  properties  of  the  ON  representation  of  E,  not  of  E  itself.  As  will  be  seen, 
it  is  desirable  to  work  with  regular  representations  of  systems  in  C.  Given  an  ODE,  it  is  a 
straightforward  matter  to  obtain  from  it  a  regular  ON  representation.  Consider  the  following 
ODE 

p1(a)w1  +  p2(cr)w2-{- ■■■pq(o)wq  =  Q  (9) 

where  pj  is  a  polynomial  in  a  —  and  Wj  are  the  manifest  variables.  Let  l  be  the  highest 
power  of  o  appearing  in  equation  (9).  Then  define  pj{cr~l)  =  pj(a)a~‘.  It  follows  that 
polynomials  pj  can  be  expressed  as 


Now  define 


:+i 

pj  =  Y,bac b>+i,i  £  Ofor  some  J 

i= 1 


A  = 


0  0  0  •••  0 

1  0  0  •••  0 

01  0  •••  0 


<E  R?x' 


0  0  •••  1  0 


B  =  [bij]  £  R'x?,  i  =  1,2,...,/ 
C  =  [0  0  -  •  - 1],  G  Rlxi 
D  =  [bt+1J]  e  Rlx? 


(10) 


(11) 


Then  the  set  of  all  solutions  w  of 


x  =  Ax  +  Bw 
0  =  Cx  +  Dw 


(12) 


is  equivalent  to  the  set  of  all  solutions  of  equation  (9).  Furthermore,  {A,  B,C,  D}  is  regular, 
by  construction.  It  is  also  minimal  in  the  sense  that  matrices  A,  B ,  C,  D  cannot  be  made 
smaller,  since  (C,  A)  is  observable  and  D  is  surjective,  a  result  proved  in  [4]. 

When  more  than  one  ODE  constitutes  a  system,  a  regular  ON  representation  may  be 
obtained  by  first  using  the  above  procedure  for  each  ODE,  and  then  using  the  interconnection 
procedure  outlined  in  the  interconnection  section  of  this  paper. 

An  output  nulling  representation  may  be  interpreted  as  an  LFT  between  a  constant  matrix 
and  an  integration  operator: 
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Figure  2:  ON  Representation  as  LFT 


where  it  is  understood  that  x  =  /„*  xdr  +  a:(0),  ie.,  the  initial  values  of  the  integrator 
states  need  not  be  0. 

The  above  LFT  interpretation  motivates  how  uncertainty  should  be  incorporated  into  our 
description,  and  leads  to  the  following  definition: 


Definition  2  An  uncertain  dynamical  system  is  in  general  a  parameterized  family  SA  = 
{SA  :  A  £  A}  of  dynamical  systems,  and  is  denoted  by  a  quadruple  EA  =  (T,W,  A,#a). 
The  parameter  A  £  A  incorporates  the  system  uncertainty. 

We  will  consider  uncertain  dynamical  systems  for  which  the  parameterized  behavior  may 


be  expressed  as 


Ba  =  {w£  L12oc  |  S(A,M)w  =  0} 


(13) 


where  M  is  a  constant  matrix,  referred  to  as  the  representation  matrix,  and  A  varies  in  a  set 
of  structured  operators,  A.  This  representation  is  denoted  (A,  M),  and  is  called  a  Generalized 
Output  Nulling  (GON)  representation  ofY*A. 


Setting  A  =  f  I  yields  the  same  behavior  as  ON  representation  {A,  B,C,  D},  where  the 
representation  matrix  M  is  partitioned  accordingly.  In  the  presence  of  uncertainty,  A  is  of 
the  form 

A  =  {diag[fl,  Au]}  (14) 

This  choice  of  A  induces  a  natural  partition  of  representation  matrix  M : 


M  = 


'  A  A  Ft 
Aun  Au  13 u 
Cn  Cu  D 


where  subscript  n  implies  nominal  and  subscript  u  implies  uncertain. 

As  in  the  case  of  no  uncertainty,  a  GON  representation  is  said  to  be  regular  if  D  is 
surjective.  We  will  be  interested  in  uncertain  dynamical  systems  for  which  there  exists  a 
regular  GON  representation.  There  exist  systems  for  which  the  regularity  assumption  is  not 
valid.  Consider  the  following  equation: 

bx  =  0  (16) 

where  6  is  a  real  parameter,  and  x  a  variable.  It  can  be  easily  verified  that  there  does  not 
exist  a  GON  representation  of  the  above  with  D  non-zero.  In  a  sense,  the  above  is  a  result 
of  choosing  a  “bad”  value  for  the  nominal  coefficient  multiplying  x,  in  this  case  0.  The 
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assumption  of  regularity  implies  that  the  system  behavior  at  A„  =  0  is  in  some  sense  similar 
to  the  system  behavior  in  a  neighborhood  of  Au  =  0.  In  the  simple  example  above,  this 
problem  could  be  remedied  by  defining  a  new  real  parameter  S  =  $  —  1,  and  rewriting  the 
above  equation  as  (1  +  S)x  =  0.  It  should  be  noted,  however,  that  the  nominal  constraint  is 
now  x  =  0,  as  opposed  to  no  constraint  on  x  when  6  =  0. 

There  are  some  transformations  on  the  representation  matrix  which  are  behavior  preserv¬ 
ing.  We  have  the  following  result: 


Lemma  1 

matrix  if 


Given  representation  matrix  M 


A  B 
C  D 


M 


T-\A  +  LC)T  T-\B  +  LD) 
PCT  PD 


,  M  is  an  equivalent  representation 

(17) 


where  L  is  any  matrix,  P  is  any  square,  invertible  matrix,  and  T  is  any  square,  invertible 
matrix  such  that  T~1AT  —  A  VA  £  A. 


The  proof  of  the  above  is  a  direct  extension  of  Theorem  3.20  in  [4].  Matrices  L,P,  and  T 
can  be  interpreted  as  follows:  L  injects  the  output  of  the  LFT  into  the  state  equation.  Since 
the  output  is  required  to  be  zero,  the  behavior  is  obviously  preserved.  Similarly,  bijectively 
mapping  the  output  of  the  LFT  by  P  to  another  output  space  does  not  alter  the  behavior, 
since  the  output  is  zero.  Finally,  T  can  be  interpreted  as  a  state  transformation. 

In  the  absence  of  uncertainty,  and  when  the  representation  is  minimal  in  the  sense  that 
matrix  M  is  as  small  as  possible,  the  implication  in  the  above  lemma  goes  the  other  way  also; 
all  equivalent  representations  are  related  as  in  equation  (17).  Whether  this  is  also  true  for 
minimal  representations  of  uncertain  systems  is  an  open  research  question,  aspects  of  which 
are  considered  in  [8]. 


4  Interconnection 

Interconnection  of  systems  without  uncertainty  in  a  behavioral  setting  is  formally  defined  in 
[2].  For  the  purpose  of  this  paper,  we  will  use  the  following  simplified  definition: 

Definition  3  The  interconnection  of  two  systems  with  the  same  signal  space  Si  =  (R+ ,  R? ,  B\) 
and  S2  =  (R+,R?,B2)  is  denoted  Si  A  S2,  where  Sx  A  S2  =  (R+,R3,j6i  0  B2) 

Thus  interconnection  can  be  interpreted  as  the  intersection  of  behaviors,  or  as  combining 
constraint  equations.  The  requirement  that  the  signal  spaces  of  two  systems  be  equal  before 
interconnection  can  easily  be  accommodated  by  augmenting  the  signal  space  of  each  system 
to  include  the  other  system’s  variables. 

The  machinery  of  interconnection  presented  below  may  be  used  when  building  a  model 
from  first  principles,  or  when  systems  are  connected  together.  Mathematically  they  are  equiv¬ 
alent,  since  both  consist  of  combining  constraint  equations.  As  stated  in  the  previous  section, 
a  system  of  ODEs  can  be  represented  by  a  regular  ON  representation  by  simply  building  ON 
representations  of  each  ODE,  and  then  interconnecting  them. 


Proposition  1  Given  regular  ON  representations  {Ax,  Bi,Cx,  Dx}  and  {A2,  B2,  C2,  D2}  of 
two  systems  Sj  and  £2  with  the  same  signal  space,  an  ON  representation  of  their  intercon¬ 
nection  Sj  A  £2  is  ,  r 


S  = 


0 


0 

a2 


B1 

B2 


Cx  0 

0  c2 


Dx 

D2 


(18) 


The  above  follows  immediately  from  the  definition  of  interconnection,  and  can  be  trivially 
extended  to  the  interconnection  of  more  than  two  systems. 

There  are  three  distinct  cases  which  may  arise.  S  may  be  dependent,  singular,  or  regular. 
A  dependent  representation  may  be  obtained  when  there  are  redundant  constraint  equations. 
In  electrical  circuits,  this  can  easily  occur  when  there  are  too  many  loop  or  node  equations.  It 
is  desirable  in  these  cases  to  simply  omit  the  redundant  equations.  Singular  interconnections 
arise  when  the  interconnection  forces  an  algebraic  constraint  between  states  in  the  subsystems. 

We  can  think  of  interconnection  in  two  contexts.  First,  when  interconnection  is  simply  an 
artifice  of  our  modeling  process,  where  we  have  broken  the  system  into  subsystems.  In  this 
case,  a  singular  interconnection  is  simply  a  flag  that  our  states  are  constrained  and  therefore 
we  might  want  to  simplify  the  model.  A  simple  example  is  the  modeling  of  two  capacitors 
in  parallel  using  two  states  for  the  voltages  across  the  capacitors.  There  exists  an  algebraic 
constraint  on  the  states  that  requires  them  to  be  equal  for  all  time. 

Second,  when  a  physical  interconnection  is  established  at  a  particular  instant  in  time.  In 
this  case,  a  singular  interconnection  would  require  that  the  states  be  matched  in  advance, 
otherwise  we  will  have  a  transient  phenomenon  which  is  not  modeled  and  is  potentially  dam¬ 
aging.  If  two  capacitors  are  connected  together  in  parallel,  one  must  be  aware  that  unless  the 
residual  charges  on  the  two  capacitors  are  compatible,  a  large  current  will  flow  at  the  time  of 
contact.  We  will  not  deal  here  with  describing  such  phenomena,  which  are  considered  in  [9]. 

We  can  extend  Definition  3  to  include  uncertainty: 


Definition  4  The  interconnection  of  two  systems  with  the  same  signal  space  x  =  (R+,  R?,  A1;  £>Al) 
and  EA2  =  (R+,R?,  A2,#a2)  is  denoted  £aiA£a2>  where  £aiA£a2  =  (R+,R?,(Ai, 

Ba3). 

A  similar  extension  of  Proposition  1  yields 


Proposition  2  Given  regular  GON  representations  (AX,MX)  and  (A 2,M2)  of  two  systems 
£ai  and  £A2  with  the  same  signal  space,  a  GON  representation  of  their  interconnection 
£ai  A  £A2  is 


Ain  0  Alnu  0  B  in 

0  A2a  0  A2nu  B2n 

J ^  0  Aln  0  Blu 

0  A-2  0  A-2  B2 

Cl.  0  Clu  0  Dx 

0  C2„  0  C2u  d2 

A  =  {diag  [//,  AUl,  AUJ  } 


(19) 
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where 


A%„  A{nu  Bin 
Mi  =  Aiun  Aiu 

.  Q  n  Ciu  I),  _ 

A;  =  {diag  [//,  AU1]  } 

In  order  to  check  whether  the  interconnection  of  two  systems  is  regular  or  not,  we  must 
start  out  with  regular  representations.  After  interconnection,  we  would  also  like  to  make 
the  resulting  representation  regular.  What  follows  is  one  step  of  a  recursive  procedure  which 
enables  us  to  construct  a  regular  GON  representation  from  one  that  is  not. 

Consider  a  GON  representation  (A ,M)  of  SA,  where  M  is  partitioned  as  in  (15).  We 
then  have  the  following  equations: 

x  =  Anx  +  AnuAuz  +  Bnw  (20) 

2  =  Aunx  +  AuAuz  +  Buw 
0  =  Cnx  +  CuAuz  +  Dw 


If  D  is  not  surjective,  there  exists  a  non-singular  matrix  P  such  that  the  last  row  of  PD 
is  zero.  Furthermore,  since  P  is  non-singular,  the  solutions  to  the  above  equations  are  not 
changed  if  the  last  equation  of  (20)  is  pre-multiplied  by  P.  Therefore  it  can  be  assumed, 
without  loss  of  generality,  that  the  last  row  of  D  is  zero,  which  results  in  the  equation 

CUrx  +  CUrAuz  —  0  (21) 

where  Cnr  is  the  last  row  of  matrix  Cn,  and  CUr  is  the  last  row  of  matrix  Cu. 

If  CHr  and  CUr  are  zero,  the  GON  representation  is  said  to  be  dependent,  which  generalizes 
our  previous  definition  of  dependent.  In  the  case  of  a  dependent  representation,  another 
representation  exists  with  at  least  one  fewer  constraint  equations.  One  such  representation 
can  obtained  by  truncating  the  last  rows  of  Cn,  Cu,  and  D.  If  the  truncated  D  is  still  not 
surjective,  the  process  may  be  reiterated. 

When  CUr  is  zero  and  Cnr  is  not  equal  to  zero,  there  exists  an  invertible  matrix  T  such 
that  CnT  =  [10  -  -  -  0] .  Letting  x  —  T~lx  results  in  the  following  equivalent  representation 
matrix  r  ,  .  ..  r  m  a 


‘  A„ 

A 

Bn' 

' T-'AnT  T~xAnu 

T-]Bn  ' 

un 

Au 

Bu 

= 

AunT 

Au 

Bu 

Cu 

D 

CnT 

Cu 

D 

We  then  have  the  following: 

x  —  Anx  -(-  Anu Au z  T  Bnw 
z  —  Aunx  +  AuAuz  T  Buw 
0  =  Cnix  +  CUiAuz  +  DiW  i  =  1,2,- 

0  =  xi 

Differentiating  the  last  equation  yields 

0  —  An,x  -(-  AnUl  Auz  T  Bniw 
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where  Ani  is  the  first  row  of  An,  AnUl  is  the  first  row  of  Anu ,  and  Bni  is  the  first  row  of  Bn. 
Note  that  the  above  equation  implies  that  aq  =  constant,  not  0,  therefore  aq  =  0  must  be 
incorporated  into  our  state  space  equations.  Therefore  representation 


M  = 


mu  Mp 

M21  Ad 


A„ 

A 

B„ 

A 

f’-un 

Au 

O 

JJU 

cni 

cUl 

P>x 

Cn 

Cu^ 

Dr— i 

Ani 

A 

finu  i 

Bni 

(25) 


along  with  the  constraint  aq  =  0  yields  the  same  set  of  solutions  w  as  M.  But  the  constraint 
aq  =  0  can  be  imposed  directly  on  representation  matrix  M  by  removing  the  first  row  and 
first  column  of  M,  resulting  in  representation  matrix  M  which  has  one  less  integrator  state. 

It  immediately  follows  that  the  above  step  could  be  repeated  until  a  regular  GON  repre¬ 
sentation  is  obtained: 


Proposition  3  Assuming  that  CUr  —  0  whenever  Dr  =  0  in  the  above  step  of  the  recursive 
reduction  procedure,  a  regular  GON  representation  can  be  obtained  from  a  non-regular  GON 
representation  M  in  at  most  N  steps ,  where  N  is  the  dimension  of  matrix  An 

If  Cu r  is  not  equal  to  zero  when  Dr  is  zero,  the  reduction  procedure  as  outlined  will 
not  work,  even  though  a  regular  representation  may  exist.  All  the  examples  which  we  have 
devised  that  cause  the  reduction  procedure  to  fail  with  a  constant  P  matrix,  however,  have 
had  no  physical  interpretation. 

The  shortcoming  of  the  procedure  is  that  multiplying  the  last  row  of  (20)  by  an  invertible, 
constant,  matrix  P  is  not  general  enough.  In  fact,  P  can  be  any  invertible  operator,  not 
necessarily  constant.  For  the  case  where  the  uncertainty  structure  consists  exclusively  of 
scalar  times  identity  blocks,  which  arise  naturally  when  dealing  with  parametric  uncertainty, 
we  have  devised  a  method  which  enables  us  to  continue  the  reduction  procedure  when  a 
regular  representation  exists,  at  the  cost  of  greatly  increasing  the  complexity  of  the  uncertainty 
structure.  The  final  representation  is  also  only  valid  locally  about  A„  =  0.  The  main  idea  is 
to  close  the  LFT  loop  on  the  uncertainty  structure,  and  reduce  D(AU)  to  row  reduced  echelon 
form  by  left  multiplying  with  matrix  P(AU),  where  Au  is  an  augmented  version  of  Au.  This 

is  t^5JV^ieiflh^<Mbwing  P  to  be  an  invertible  operator,  it  still  may  not  be  possible  to  make 
CUr  zero  when  Dr  is  zero.  For  parametric  uncertainty,  this  occurs  when  D(AU)  loses  rank 
at  Au  =  0,  but  not  in  a  neighborhood  of  Au  =  0.  An  interpretation  is  that  there  may 
exist  algebraic  constraints  on  the  integration  states  at  Au  =  0,  but  not  in  a  neighborhood  of 

Au  =  0.  A  simple  example  is  depicted  in  Figure  4. 

The  nominal  resistance  is  zero,  which  implies  that  the  voltages  across  the  capacitors  must 

be  nominally  equal,  but  may  be  different  for  all  other  values  of  Sr.  Clearly  the  problem  arises 
due  to  the  zero  nominal  resistance,  a  poor  modeling  choice. 
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Figure  3:  Multiplying  by  P(AU) 


Vi  o 


o  P2 


Figure  4:  System  without  Regular  Representation 


Issues  such  as  whether  physically  motivated  examples  exist  for  which  the  procedure  breaks 
down,  how  to  deal  with  full  uncertainty  blocks  when  using  the  augmented  reduction  procedure, 
and  how  to  reduce  the  complexity  of  the  augmented  procedure  are  open  research  questions. 

Non-regular  representations  which  can  be  reduced  to  regular  representations  will  be  called 
singular,  which  generalizes  the  previous  definition  of  singular. 

5  Input-Output  Map  Representation 

As  discussed  in  [1],  it  may  be  required  to  represent  the  behavior  of  a  system  SA  in  input- 
output  form.  This  may  be  the  case  if  existing  robust  control  methodologies  are  to  be  used  for 
the  purpose  of  designing  robust  controllers,  as  these  methods  require  a  distinction  between 
system  inputs  and  outputs.  The  reader  is  referred  to  [2]  for  a  characterization  of  input  and 
output  variables.  We  will  restrict  ourselves  to  input-output  maps  which  are  nominally  proper. 

Given  a  GON  representation,  we  need  to  know  how  many  variables  may  be  chosen  as 
inputs. 

Proposition  4  Given  a  regular  representation  (A,  M)  with  D  £  Rrx?,  any  input-output  map 
must  have  r  outputs  and  q  —  r  inputs. 

Proof: 

First,  note  that  the  input-output  map  must  be  valid  for  the  nominal  system.  Thus  we 
can  restrict  ourselves  to  the  case  of  no  uncertainty.  Then  it  is  a  simple  matter  to  convert  the 
ON  representation  to  an  AR  representation,  and  then  invoking  Theorem  VIII.7  in  [2]  (which 
states  that  an  AR  representation  with  r  independent  equations  must  have  r  outputs). 


10 


□ 


Let  (A,  M )  be  a  regular  representation  of  SA-  We  then  have  the  following: 

p  =  AAp+Bw  (26) 

0  =  CAp  +  Dw 


Since  D  is  surjective,  there  exists  permutation  matrix  II  =  [n1n2]  such  that  HIl!  is  invertible. 


Let  w  =  II 


y 

u 


Then  we  have 


p  =  AAp  +  Bxy  +  B2u  (27) 

0  =  CAp+D1yJrD2u 


where 


bx  =  Bnlt  b2  =  b n2,  Dy  =  d nlt  d2  =  d n2 


Isolating  y  and  substituting  in  the  p  equation  yields 

p  =  (A  —  ByDi~1C)Ap  + (B2  —  BiDi~1D2)u 
y  =  -Dy~lCAp  -  Dx~xD2u 


(28) 


Note  that  the  above  is  in  the  form  of  an  LFT  relating  y  and  u.  As  of  yet,  all  operations 
performed  have  been  behavior  preserving.  If,  however,  it  is  required  that  u  vary  freely  and 
determine  y ,  it  must  be  possible  to  solve  for  p,  ie.,  for  a  given  u  £  Ll?c,  there  must  exist  a 
p  £  Ll2c ■  This  might  not  be  the  case  for  all  possible  values  of  A  in  A.  We  have  the  following 
main  result: 


Theorem  1  Given  regular  GON  representation  (A,  M)  ofT, A  =  (R+,  R9,  A,  B&),  there  exists 
a  partition  of  R?  =  Rp  X  Rm,  a  permutation  matrix  II,  a  constant  matrix  Mj/o,  and,  a  number 
e  strictly  greater  than  zero  such  that 


=  {n 


y 

u 


y  £  L‘2°°(RP),  u  e  40C( Rm)  I V  =  ^(a,  miio)u } 


(29) 


for  all  A  £  B A  e .  (A ,M//o,H)  will  be  referred  to  as  an  I/O  representation  of  ■ 


Proof:  Consider  the  following  equations 

x  —  Anx  T  AnuAuz  +  Bnu  (^0) 

z  =  Aunx  +  AuAuz  +  Buu 
y  —  Gri x  Gu  Auz  T  Bj/qu 

which  correspond  to  (28)  when  the  matrices  are  partitioned  into  their  corresponding  nominal 
and  uncertain  parts,  and  p  —  (x,z).  It  is  clear  that  for  sufficiently  small  Au,  (/  -  Au A„)  will 
be  invertible.  In  other  words,  there  exists  an  e  >  0  such  that  (/  -  Au Au)  is  invertible  for 
||  Au  || <  e.  We  then  have  the  following 
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where 


x  =  Lx  +  Qx  +  / 


(31) 


L  =  An 

0  =  AnuAu{I  ~  AUAU)~' Aun 
f  =  (AnuAu(I  —  Au  Au)  +  -Bn)^ 

We  have  the  following  technical  lemma,  whose  proof  may  be  found  at  the  end  of  this  section, 
which  is  central  to  proving  the  theorem: 

Lemma  2  Given  f  G  L2oc,  0  :  Llfc  —*  Llfc,  ||  0  ||<  M  <  oo,  0  causal  and  linear,  and  L  a 
constant  square  matrix,  V  initial  conditions  x(0)  =  x0  3  unique  x(t)  for  t  >  0  such  that  x{t ) 
solves  (31 )  a.e.. 

Invoking  Lemma  2,  we  conclude  that  x  can  be  determined  from  u,  and  that  x  G  L!fc.  It 
then  follows  from  (30)  that  z  can  also  be  solved  for,  and  that  2  G  Llfc.  Finally,  since  y  is  a 
bounded  function  of  x,  z,  and  u,  y  must  also  be  in  Lir- 

□ 


We  then  have  the  following  corollary: 


Corollary  1  A  valid  I/O  representation  matrix  Mj/0  of  SA  is 


Mjjo  = 


Ai/o  Bj/o 

0 

b 

0 

1 _ 

A  -  B1D1~1C 
-Dl~1C 


b2-b1d1~1d2 

-D{~lD2 


(32) 


How  can  we  get  a  good  bound  for  e?  One  way  is  to  perform  a  p  test  on  the  matrix  Au 
with  block  structure  Au.  By  the  definition  of  p,  see  [10],  e  =  ■ — . ~t . -r  is  the  largest  value  for 

/I  A  u  V-™  u  ) 

which  (/  -  Au  A„)  is  invertible  for  ||  Au  ||<  e.  We  then  have  the  following  corollary: 


Corollary  2  A  value  of  e  consistent  with  Theorem  1  is  ^  ^  )f  where 


Aijo 


A  A 


(33) 


Proof  of  Lemma  2:  Pick  an  interval  [0,  T],  where  T  is  to  be  chosen.  Define  x0(t)  =  x0,  and 

£n+i(t)  =  z“o+  f  {Lx„(r)  +  (0ain)(r)  +  /(r))dr  ,  n  >  0,  0  <  t  <  T  (34) 

J  0 


For  n  >  1 


Vn+i(t)  ~  Xn(t)  \  <  f  (|  L(xn(r)  -  £„_i(t))  |  +  |  (Q(xn  -  s„_i))(r)  \)dr 
Jo 

—  d(L)t  2  II  xn  —  Xn-l  ||2  [oit]  +  t  2  ||  ®{xn  —  xn- 1)  ||2  [o,t] 

<  {a(L)  +  M)ti  ||2[M  (35) 


which  implies 
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(36) 


T 

||  xn+ 1  —  xn  II 2 [0,T]  —  ^)ll  Xn  ~  xn-l  ll2[0,T] 

Now  choose  T  to  make  ^(d(L)  +  M)  <  1.  Then  consider 

OO 

x(t)  =  X0(t)  +  -  xn(t))  6  X[°’T]  (37) 

n=0 

Then  it  follows  that  1o,t] 

xn(t)  -2-  x(t)  (38) 

Taking  limits  on  equation  (34),  and  using  the  continuity  of  0  and  integration  operators  in 

4°'T],  yields  t 

x(t)  =  x0+  I  (Lx(t)  +  (Qx)(t)  +  /(r))dr  (a.e.)  (39) 

Jo 

which  implies  that  x  solves  (31)  a.e.  in  [0,T].  The  argument  above  can  be  easily  extended  to 
show  uniqueness. 

It  remains  to  see  whether  we  can  extend  the  solution  x  to  the  real  line.  Clearly,  it  suffices 
to  show  we  can  do  so  in  the  interval  [0,2T],  Define 


(0  ,T  <  i  <  2T 

(40) 

and 

Xn+1  (t)  =  Xo  +  /  (Lx„(t)  +  (exn){r)  +  f(r))dr  ,  n  >  0,  0  <  t  <  2T 

J  o 

(41) 

Now, 
(41)  in 

for  t  <  T,  a;n+i(t)  =  x(t).  This  uses  the  fact  that 

[0,T],  and  that  0  is  causal.  Then  we  have,  for  T  <t<2T  and  n  >  1, 

x(t)  solves 

I  ®»+i (0  -  Xn(t)  |<  (d(X)  +  M) ||  xn  -  xn_i  ||2[Ti<](/  -  T)i 

(42) 

At  this  point,  the  same  argument  as  before  applies,  and  we  have  a  solution  x{t)  on  [0,2T]. 

□ 


6  Example 

The  following  example  illustrates  some  of  the  various  techniques  and  algorithms  presented 
in  the  previous  sections.  It  consists  of  the  interconnection  of  two  mechanical  systems  with 
uncertain  parameters.  Define  the  parameters  as 

&i  =  1  +  &kn  bi  =  1  +  Stl,  m  i  =  1  +  6mi 

k-2  =  1  +  Sk2,  b2  =  0  +  «562,  m2  =  2  +  8m2 

The  uncertainty  can  take  on  various  forms.  For  example,  the  damping  of  system  1  could 
be  interpreted  as  being  nominally  1,  with  a  constant  or  time  varying  perturbation  to  the 
nominal  damping.  It  could  also  be  interpreted  as  being  nominally  1  with  the  addition  of  an 
operator  that  acts  on  the  velocity  of  the  mass  and  converts  it  to  a  force. 

The  equations  of  motion  governing  the  two  systems  are 

/i  =  (1  +  8mi)xi  +  (1  +  #6i)*'i  +  (1  +  ^*i)®i  (43) 

/2  =  (2  +  Sm2)x2  +  8b2xi  +  (1  +  8k2)x2 
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/ 

/ 


h 

— wv- 

!  xi 

- ! 

1 

i 

>- 

A 

— - 

lTi\ 

/  /  7 


Z 


Z  WZ 


SYSTEM  1  Figure  5:  Systems  1  and  SYSTEM  2 


In  order  to  construct  GON  representations  of  the  two  systems,  we  need  to  modify  the  proce¬ 
dure  outlined  in  the  representation  section  to  allow  for  uncertainty.  A  GON  representation 
of  system  1  can  be  constructed  as  follows: 

Define 

Zi  =  xx  (44) 

z2  =  Xi 
Zz  =  Sklzi 
Z4  —  SblZ2 

z5  =  (1  +  ^mi)  fimi(Zl  +  Z2  +  Z3  Z4  —  fi) 

=  firm  ( Z\  -(-  Z2  +  z3  +  z4  —  Z5  —  /i) 


then 


z2  —  Z5  —  (z  j  +  Z2  +  Z3  +  Z4  —  /i) 


We  then  have  the  following  equations 


where 


p  =  AxAxp+  B1wl 
0  —  C\  Axp  -(-  D\W\ 


z  =  [z1,z2,z3,z4,z5]  =  AjP 


Ai 


C\ 

Ax 


w  i 


'  0 

1 

0 

0 

0 

'0 

0  ' 

-l 

-1 

-1 

-1 

1 

0 

1 

l 

0 

0 

0 

0 

,BX  = 

0 

0 

0 

1 

0 

0 
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0 

0 

1 

1 

1 

1 

-1 

_0 

-1  _ 

1  0 

0  0 

0 

,DX 

=i 

-1  o] 

/ 


hrfik1,fib1,fin 


]} 


[*l>/l] 


(45) 

(46) 


(47) 
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which  gives  us  the  required  regular  GON  representation  of  system  1.  A  similar  construction 
can  be  performed  for  system  2,  yielding 


‘010  0  o' 

o 

o 

_ 1 

=1  0  =±  1 

0  1 

1  0  0  0  0 

,B2  = 

0  0 

0  10  0  0 

1  0  1  i 

L  4  4  4  2  j 

o  o 

Hi  ° 

1 _ 

c2  =  [ioooo],i>2=[-io] 

A  2  =  \^diag 

w  =  [*2,/2] 


/ 


/2,  Vm2 


(48) 


Let  us  define  the  interconnection  of  the  two  systems  as  follows: 


/ 


Xi  =  x2 


/ — 'Vvv 

/— irO- 

/  - - - 7  y— 


mi 

m2 

I — vAAA 

Qt: 


\ 

K 


Figure  6:  Interconnected  system 


*i  =  a:2 

/i  +  /2  =  / 


(49) 


where  (49)  has  a  trivial  regular  GON  representation.  This  corresponds  to  “bolting”  masses 
1  and  2  together.  Combining  systems  1  and  2  with  the  interconnection  constraints  results  in 
the  following  GON  representation  of  the  interconnected  system: 
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t  >  ^fc2  5 


2 »  “m2 


(50) 


Note  that  the  interconnection  is  not  regular,  since  D  is  not  surjective.  This  is  expected, 
since  we  are  forcing  the  two  masses  to  move  together,  and  thus  their  positions  and  velocities 
are  constrained  to  be  equal.  Thus  we  would  expect  the  system  to  have  only  two  integration 
states,  not  four. 


After  performing  the  reduction  procedure  (which  terminates  after  two  iterations),  we  have 
the  following  regular  GON  representation  of  the  interconnected  system: 
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Using  Proposition  4  and  Theorem  1,  this  system  must  have  one  input,  which  is  restricted  to 
be  one  of  /,  /j,  or  f2.  Choosing  /  as  an  input,  Corollary  1  implies 
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Several  things  about  the  I/O  map  are  worth  noting.  First  of  all,  the  eigenvalues  of  A//on 
are  —.1667  ±  .7993?.  These  correspond  to  the  natural  modes  of  the  unperturbed  system,  as 
can  be  easily  verified.  Secondly,  if  we  assume  that  the  perturbations  are  real  constants,  a  Real 
H  test  on  the  uncertainty  portion  A//ou  of  A//o  yields  a  result  of  .6667.  This  indicates  that 
there  exist  perturbations  of  size  1.5  which  will  cause  the  I/O  representation  to  be  ill-posed. 
This  corresponds  to  Smi  and  Sm2  being  —1.5,  which  would  make  the  effective  mass  of  the 
interconnected  system,  irij  +  m2,  equal  to  0.  It  is  interesting  to  note  that  a  Real  fi  test  on 
system  1  before  interconnection  would  have  yielded  a  value  of  1.  This  of  course  corresponds 
to  Smi  being  —1,  and  causing  the  mass  to  be  0.  The  interconnection  thus  extends  the  region 
of  validity  of  the  mass  uncertainty.  This  might  flag  that  something  non-physical  might  occur 
in  the  interconnected  system  for  a  valid  value  of  the  uncertainty.  For  the  above  example,  this 
reflects  the  fact  that  we  don’t  expect  mass  1  to  be  negative,  even  though  the  interconnected 
system  is  still  well  posed  for  a  negative  value  for  mass  1. 


7  Conclusions 

In  this  paper,  an  attempt  has  been  made  to  relate  robust  control  and  behavioral  paradigms 
by  incorporating  structured  uncertainty  into  the  description  of  behavioral  systems.  This  is 
accomplished  by  expressing  behavioral  systems  as  kernels  of  operators  obtained  by  an  LFT 
between  a  constant  matrix  and  an  uncertainty  structure.  An  uncertain  dynamical  system 
has  been  defined,  and  methods  of  interconnection  developed.  Finally,  a  method  for  obtaining 
input-output  maps  from  behavioral  descriptions  of  uncertain  dynamical  systems  has  been 
outlined,  which  makes  possible  the  use  of  existing  robust  control  methodologies. 
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